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1. INTRODUCTION 

Nowadays, economic-emission problem is getting considered in today’s electricity market and it is 
very crucial to control the emission resulted. The good quality power system will bring the great satisfaction 
to the consumer since most of the daily routine depends on electricity. Economic Load Dispatch (ELD) has 
become the critical problem in electrical power generation systems. The aim of ELD is to find an optimal 
power so that the generation cost will minimized while satisfying load demand and all other operational 
constraints [1]. 

Operating at minimum cost can no longer be considered alone since the consumer need a clean 
environment [2]. The environmental aspect need to be taken into consideration when burned the fossil fuels 
at thermal power plant. The emission released such as nitrogen oxide (NOx), carbon dioxide (CO2), Sulphur 
dioxide (SO2) and etc. are hazardous to the environment. This emission has become the main cause of air 
pollution [3]. The US Air Act Amendment of 1990 mandates that the electric utility industry should reduce 
its SO2 emissions by 10million tons/year and the NOx emission by 2 million tons/year from the 1980 level in 
order to reduce the environmental pollution [4]. It shows that the use of fossil fuels to generate electricity 
need to be minimized and controlled due to environmental pollution of burning fossil fuels [5, 6-8]. 
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There are many different methods have been used as a mathematical tool for optimization 
problems. It can be divided into two categories which are traditional and metaheuristic methods. The 
Metaheuristic methods are recently used for many applications since it capable to find global optimal 
solution for economic-emission load dispatch (EELD) problems [1, 9-19]. In this project, Firefly Algorithm 
(FA) is applied to investigate the optimal cost and emission, robustness and convergence characteristics of 
EELD problem. The FA is first introduced by Xin-She Yang in 2008 and the behavior of the algorithm is 
inspired by fireflies’ movement [3]. Moreover, the main advantages of FA are simple concept and easy to 
implement since it is based on the global communication among the swarming particles [20-22]. 


2. PROBLEM FORMULATION 

Sustaining load demand while minimizing emission level and total generation cost requires big 
challenge as a multi-objective optimization problem. The elements of the load dispatch problem that include 
the objective functions and operation constraints are deeply explained as follows: 


2.1. Economic Load Dispatch 
The aim of Economic load dispatch (ELD) is to reduce the total generator cost while sustaining load 
demand and operational constraints. The formulation of ELD problem defined as [23]: 


ng ng 
Fo =F (R)= 1a, +bP +¢,P? 
i=l i=l 


() 


where, Fi (Pi) is generation cost function ($/h), ai, bi, ci are the cost coefficients for the ith 
generator, ng is the total number of dispatchable plants and Pi is the power output of ith plant. 

The objective function of ELD with valve point effect can be represented as the sum of a quadratic 
and a sinusoidal function as follows [23]: 


ng ng : 
Fo => F(P)= > (4, +52 +¢,P’)+ le, sin(f(P™ -P))| 
i=l i=l (2) 


where, ei, fi are the cost coefficients for ith generator considering the valve-point effect. 


2.2. Emission Load Dispatch 
Emission load dispatch is to reduce the total stack emission of nitrogen oxide (NOx), sulphur oxide 
(SOx) and carbon dioxide (CO2) for the entire system. The objective function defined as [19]: 


ng ng 
Fy => E(P) =a + BP +Y, 
ial ial (3) 


While, the objective function of emission load dispatch considering valve-point effect can be 
represented as follows [19]: 


ng ng 
F, = DE (P)=)1a,P? + BP +7 + (nxexp(6 x P)) 
i=l i=l (4) 





a,, 


where, Ei (Pi) is the total emission released in kg/h, i?Y7 are the emission coefficient of ith 


generating unit and 7,5; are the emission coefficients of ith generator with valve-point effect. 


2.3. Combined Economic and Emission Load Dispatch: Weighted Sum Method 

The aim of combined economic and emission load dispatch to minimize total generator cost and 
emission level that produce energy in thermal power plant simultaneously. To handle this problem, the 
weighted sum method is used to find the set of optimal solutions between these two conflicting objectives 
based on weight factor (w) as follows [5]: 
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F= Ww) K(P)+(1—wyh, E(B) 
i=l i=l (5) 


The EELD problem is expressed in a single objective form that are formulated by combining two 
independent objectives which are total cost and emission level. The weighting factor (w) can be any number 
between 0 and | based our references. The price penalty factor hi is the ratio between maximum cost and 
maximum emission level of corresponding generator in $/kg [4, 19]: 


rf 7 F.. (P™*) 


EG) 4= 1,2, cue (6) 


There following steps can be used to determine the price penalty factor for a particular load 
demand [19]: 

Calculate the ratio between maximum generation cost and maximum emission level “hi” of 
each generator. 

The values of “hi” must arrange in increasing value order. 

Add the maximum generated power of each generator one by one starting with the smallest unit 


Peano 
price penalty factor until Le ea ‘ 


At the point, hi associated with the final unit in this process is the nearest price penalty factor 
equivalent for the given load. 


2.4. Problem Constraints 

There are two operational constraints will be used in this problem constraints such as equality and 
inequality. The power balance constraint is referred to an equality constraint which expressed the total 
generated power must equal the total power demand and total power losses in the system. The equality 
constraint is formulated as follows [24]: 


= (7) 


where, PL is the total real power losses and PD is the total power demand. 
min 
The generating power limits is the inequality constraints based on minimum ( F ) and maximum 


("' ) power output as follows: 


PE apap 
i I i (8) 

The transmission network system losses are represented as function of the real power and the B- 
coefficient matrix using Kron’s loss formula as follows [24, 19]: 


ng 


ng ng 
P, = De bs (PB,;P,) - » BoP, + Boo 
i=l i=l zy 


(9) 


where, Bij is the loss coefficient matrix. 


3. FIREFLY ALGORITHM FOR SOLVING ECONOMIC-EMISSION LOAD DISPATCH 

The main purpose for a firefly’s glow is to act as a prompt system to attract other fireflies. This 
algorithm is formulated by assuming three idealized rules [25]. The fireflies attractiveness of the fireflies is 
related to the light intensity and defined as follows [2]: 


Goa Bie. sam (10) 


where, By is the attractiveness at r =0 and r is the distance between two fireflies. 
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The movement of Firefly Algorithm is defined as the following: [26] 


t+l t t t 
Me SKA Ah, HH Eas, 


(1) 


where, x is the position of each fireflies, Fis the attractiveness of the ith firefly receives from the jth 
t 


fireflies, “is the randomization parameter and ®i ig the vector of random numbers drawn from Gaussian 
distribution on the tth iteration. 

This section described about the implementation steps of Firefly Algorithm for combined economic 
and emission problem. The effectiveness of FA will be compared with other algorithms. The steps are 
implemented as follows: 


Step 1: Define the objective function of economic-emission dispatch problem using (1) and (3). 
Step 2: Set the parameter setting of FA which are a=0.3, /=0.1 and y=1 and number of population is set to 
30. 
Step 3: Insert randomized with respect to maximum power generator (Pax) and minimum power generator 
(P, min) 
Step 4” Evaluate the new solution of Firefly Algorithm with respect to objective function to consider losses. 
The implementation step for each objective function are: 

a. For economic dispatch, (1) is used to minimize the cost of IEEE 118-bus 14-unit test system. 

b. For emission, (3) is used to minimize the emission dispatch of IEEE 118-bus 14-unit test system. 

c. For minimization of combined economic and emission load dispatch, (5) is used. The weighting factor 

are selected from 0 to | with steps of 0.1. 

Step 5: Rank the new solutions from maximum to minimum depends on the objective function value and 
store the best solution. 
Step 6: The constraints in (7), (8) and (9) can be satisfied by using Modification of Infeasible Particle (MIP) 
approach in [27]. 
Step 7: Get the movement of firefly which get the optimal solution using (10) and (11). 
Step 8: Store the best solution respected to the maximum iterations. Step 4 is repeat until reached maximum 
number iterations. 
Step 9: The best solutions are stored respected to the number of iteration using Fuzzy-based Mechanism. 
Step 10: End of algorithm. 


3.1. Best Compromise Solution: Fuzzy-Based Mechanism 
The best compromise solution is selected among pareto optimal set can be determined using Fuzzy- 
based approach. Each of objective function is represented by a mechanism function 4; as shown in (12) [28]. 


1 F; < Be 
Finex_rp, 
Le J J Fmin <F < pmax (12) 
J fot tenis md re ad) 
0 F; > Eines 


where, pee Fi” = minimum and maximum value of the j-th objective function. 
A normalized mechanism function need to be evaluated for every solution i-th using (13), [28]. 


pep aie (13) 


The highest value of is chosen as the best compromise solution. 


4. SIMULATION RESULTS 

The objective of this result is to develop Firefly Algorithm for solving CEED problem using 
Weighted Sum Method and Fuzzy-Based Mechanism. The simulation and analysis were conducted to 
identify the performance of Firefly Algorithm on IEEE 118-Bus14-unit test system with total power demand 
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is 950 MW. The parameters of Firefly Algorithm take the following values for all test system which are 


a=0.3, B=01 and¥ =!. The maximum iteration and number of population are set to 1000 and 30, 
respectively. The generator data of cost and emission ARE taken from [28]. 


4.1. Minimizing Cost and Emission Individually of IEEE 118-Bus 14-Unit System 

To validate the effectiveness of the proposed FA, each cost and emission objective function is 
solved individually. The results obtained this test system are shown in Table 1. It highlights the minimum 
cost and emission are 4303.6009 $/h and 25.2374 kg/h respectively. Table 2 presents the comparison of the 
FA with other method. The results show that the FA has achieved significantly better results compared 
to MHSA. 


Table 1. Optimal Cost and Emission Obtained by FA 











Generators Minimize Cost Minimize Emission 
P, (MW) 104.6483 70.8384 
P, (MW) 92.7271 50.0000 
P; (MW) 50.0000 78.0791 
Py (MW) 50.0000 88.8467 
P; (MW) 50.0000 67.7914 
Ps (MW) 50.0000 50.0000 
P; (MW) 50.0000 73.4668 
Ps (MW) 50.0000 72.3696 
Py (MW) 59.4909 73.7123 
Pio (MW) 64.6067 90.0105 
Pi, (MW) 64.1038 50.0000 
Piz (MW) 176.9597 72.5723 
Pi3 (MW) 50.0000 72.3089 
Pi4 (MW) 50.0000 50.0000 
Total Power (MW) 962.5366 959.9961 
P, (MW) 12.5366 9.9961 
Pp (MW) 950.0000 950.0000 
Total Cost ($/h) 4303.6009 4547.9254 
Total Emission (kg/h) 394.1541 25.2374 





Table 2. Comparison for Minimize Cost and Emission Individually between the Methods for Test System 2 








Algorithm Minimize Cost Minimize Emission 
Total Cost Total Emission (kg/h) Total Cost Total Emission (kg/h) 
($/h) ($/h) 

MHSA [28] 4304.9500 357.3390 4539.2280 27.8920 

Proposed FA 4303.6009 394.1541 4547.9254 25.2374 





Figure | shows the convergence characteristics for the best result of cost and emission. It can be 
seen that the approximate optimal values reached at around 200 and 300 iterations for cost and emission 
minimization respectively. The figure proves the high-speed convergence of FA in solving EELD problem. 


4.2. Combined Economic and Emission Load Dispatch of IEEE 118-bus 14-Unit System 

The results obtained by FA for solving combined economic and emission dispatch are tabulated in 
Table 3 for different weight values. From the table, it clearly shows the value w=0 is minimize the emission 
while value w=1 is minimize the total generation cost. 


Table 3. Pareto front Solutions 








w Cost Emission Lc be a 

0.0 4547.9254 25.2374 0.0000 1.0000 0.0714 
0.1 4522.6029 26.5943 0.1036 0.9964 0.0786 
0.2 4495.9941 31.8206 0.2125 0.9824 0.0854 
0.3 4466.5354 43.5049 0.3331 0.9511 0.0918 
0.4 4428.8999 66.2426 0.4871 0.8901 0.0984 
0.5 4387.9832 103.7431 0.6546 0.7897 0.1032 
0.6 4352.3001 146.6544 0.8006 0.6747 0.1054 
0.7 4330.8089 186.7546 0.8885 0.5673 0.1040 
0.8 4313.6995 245.1099 0.9586 0.4110 0.0979 
0.9 4305.3779 282.4456 0.9926 0.3109 0.0932 
1.0 4303.5753 398.5150 1.0000 0.0000 0.0715 
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Based on Figure 2, the total cost to generate electricity at thermal power plant is increased as the 
weighting factor is decreased from 1 to 0. Meanwhile, for the total emission released from the generator is 
decreased when the weighting factor decreased. It is shown that the total cost and emission are conflicting to 
each other. Therefore, the fuzzy based method in section 3.1 is used to determine the best compromise 
solution among the solutions in Figure 2. 

Therefore, the optimal power output for best compromise solution obtained by fuzzy approach is at 
w=0.6. The total cost obtained for the best compromise solution is 4352.3001 $/h whereas the total emission 
produced is 146.6544 kg/h which correspond to w=0.6. The results obtained by proposed FA has been 
compared with the results of other algorithms such as NGSA-IT, MHSA and HAS as shown in Table 4. It is 
obvious that the total cost achieved by proposed FA is significantly better than NGSA-II and MHSA in terms 
of cost while the emission level is slightly lower than HSA. 


450 























= Total Cost a 
—Total Emission =z a 20 
zB 49 200 2 2 an 
5 2 =F 250 
= 4370 0 5 6 (150 
= 100 
50 
4250 0 
0 200 400 600 800 1000 0 | ; 
Number of Iterations 4300 4350 4400 4450 4500 4550 
Total Cost (S/h) 
Figure 1. Convergence Characteristics of FA for Figure 2. Pareto Front Obtained by FA with WSM 


minimizing Cost and Emission Individually 


Table 4. Comparison of the Best Compromise Solutions Obtained by Different Methods 





Algorithm NGSA-II [29] MHSA [29] HSA [29] _ Proposed FA 
P; (MW) 69.4618 79.3683 99.2872 98.9649 
P, (MW) 65.9624 62.1868 68.1241 64.1944 
P; (MW) 51.9713 50.0000 63.5609 50.0000 
Py (MW) 72.3437 61.1287 62.0231 67.3355 
P; (MW) 73.6853 78.6284 56.3413 63.5035 
Ps (MW) 70.3004 51.5320 50.0000 50.0000 
P; (MW) 67.5641 67.3322 50.0000 50.0000 
Ps (MW) 69.5526 50.0000 50.0000 50.0000 
Py (MW) 69.4421 79.3041 82.4229 85.1549 
Pio (MW) 81.1416 93.6602 100.2080 99.0745 
Pi, (MW) 71.3414 68.3105 51.5145 58.2738 
Pi. (MW) 70.7760 96.0040 126.6166 123.4539 
Pi3 (MW) 69.9194 68.9558 50.0000 50.0000 
Pi4 (MW) 68.5146 53.2431 50.0000 50.0000 
Pp (MW) 900.0000 900.0000 900.0000 900.0000 
Total Cost ($/h) 4534.3995 4432.298 4326.4785 4352.3001 
Total Emission (kg/h) 122.5034 100.708 153.7821 146.6544 





5. CONCLUSION 

In this paper, the performance of Firefly Algorithm with fuzzy approach has been investigated to 
optimize the total generation cost and emission level in power dispatch. The proposed method used weighted 
sum method to combined both objective function (cost and emission level) and generate the set of possible 
solution based on various weighted values. From this solution, fuzzy mechanism is utilized to determine the 
best compromise solution for minimizing cost and emission of power generation simultaneously. The 
effectiveness of proposed algorithm has been tested on IEEE 118 nus with 14-unit system by using 
MATLAB software. From the simulation results, it highlights the effectiveness of proposed method in terms 
of optimal cost and emission, convergence characteristics, Pareto front solution and comparison of the best 
compromise solution. It found that FA affirms its good performances to obtain better optimal cost and 
emission for CEELD problems. Therefore, it can be concluded that FA confirms the effective high-quality 
solution for solving CEELD problem and can be used for solving any single and multi-objective problems. 
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